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We study the behavior of droplets for two dimensional Ising spin glasses with Gaussian interac- 
tions. We use an exact matching algorithm which enables study of systems with linear dimension 
L up to 240, which is larger than is possible with other approaches. But the method only allows 
certain classes of droplets to be generated. We study single-bond, cross and a category of fixed 
volume droplets as well as first excitations. By comparison with similar or equivalent droplets gen- 
erated in previous works, the advantages but also the limitations of this approach are revealed. In 
particular we have studied the scaling behavior of the droplet energies and droplet sizes. In most 
cases, a crossover of the data can be observed such that for large sizes the behavior is compatible 
with the one-exponent scenario of the droplet theory. Only for the case of first excitations, no clear 
conclusion can be reached, probably because even with the matching approach the accessible system 
sizes are still too small. 



PACS numbers: 75.50.Lk, 02.60.Pn, 75.40.Mg, 75.10.Nr 



I. INTRODUCTION 

Ising spin glasses are amongst the most-frequently 
studied systems in statistical physics 1 '. However, despite 
more than two decades of intensive research, many prop- 
erties of spin glasses, especially in finite dimensions, are 
still not well understood. For two-dimensional Ising spin 
glasses it is now widely accepted that no ordered phase 
for finite temperatures exist For d = 2 the be- 
havior is usually described by a zero-temperature droplet 
scaling (DS) approac h 7 : 8 : 9 . 

The droplet picture assumes that the low-temperature 
behavior is governed by droplet-like excitations, where 
excitations of linear spatial extent I typically cost an en- 
ergy of order I 6 . These excitations are expected to be 
compact and their surface has a fractal dimension d s < d, 
where d is the space dimension. Furthermore it is usu- 
ally assumed that the scaling behavior of the energy SE 
of different types of excitations, e.g. droplets and do- 
main walls, induced by changing boundary conditions, 
are described by the same exponent 9. 

The value of 9 for domain walls in systems with 
Gaussian disorder has been determined by several 
studiesS^iSiiSiiiiiS with results close to the most recent 
result 9 = —0.287(2). For simplicity, we will assume 
9 = —0.29 from now on. On the other hand, in a pre- 
vious study of droplets by Kawashima and Aokiia using 
Monte Carlo simulations of moderate sized systems, a 
different exponent 9 — —0.45(5) was found. Note that 
droplets can be generated by using several methods, and 
the resulting scaling behavior may in principle depend 
on the droplet type considered. By using exact ground- 
state (GS) algorithms^, so called matching algorithms 
(see Sec. |Hj, one can study much larger system sizes than 
is possible with other approaches. On the other hand 



matching algorithms are less flexible than Monte Carlo 
methods, and not all types of droplets can be generated 
using them. But recently, through the application of an 
extended matching algorithmic, we were able to gener- 
ate the same type of droplets as in Ref . Q but for much 
larger systems. For large systems, the value 9 = —0.29 of 
the exponent was again found. To be precise, the scaling 
of the droplet energies exhibited a crossover of the form 



AE = Af + BV 



(1) 



with 9 = —0.29, to = 1, while for small sizes, the be- 
havior was compatible with an apparent exponent of 
9' = —0.47(1). Hence the predictions of DS have been 
verified for this type of droplet, simply by studying larger 
system sizes. 

Recently, the same crossover phenomenon was sought 
by Berthier and Young— when studying droplets, ob- 
tained by a Monte Carlo minimization in a system with 
size L = 64 while fixing the size of the droplets to either 
one value or to a range of values. In both cases, for the 
largest droplets, the scaling was determined by an expo- 
nent resembling the domain-wall exponent, i.e. close to 
9 — —0.23. But only for the second type of droplets with 
a range of sizes could a crossover phenomenon be ob- 
served, and then it only occurred for very small droplets. 
For a third type of droplet, no crossover and a scaling 
exponent 9 — —0.32 close to the domain-wall value was 
founcUC. 

On the other hand, Picco, Ritort and Salesi2i2S have 
studied the lowest excitation in systems with sizes up to 
L = 16 using a transfer matrix approach. Through a scal- 
ing argument connecting lowest excitations to droplets 
(see Sec. [Vj, their result implied 9 = —0.47. 

In order try to clarify the above mentioned results and 
relate them to our previous results^ we study in this 
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paper different types of droplets by applying again the 
matching algorithm. This allows us to go to larger sizes 
than in the studies where no matching algorithm was 
applied. On the other hand, the matching algorithm is 
less flexible and enables us to generate only certain types 
of droplets. Where differences in the scaling behavior 
arise amongst the various droplet types we shall attempt 
to understand their origin. 

The Hamiltonian which is studied here is the usual 
Ising spin glass model: 



H 



(2) 



where the sites i lie on the sites of a square lattice with 
N = L 2 sites, Si = ±1, the Jy have a Gaussian distribu- 
tion of zero mean and unit variance and couple nearest- 
neighbor sites on the lattice. 

In the next section, we will outline how the matching 
algorithm works and how the droplets are generated. In 
the main four sections of this paper, four different types 
of droplets are studied, their scaling behavior is investi- 
gated and compared to previous results, if available. In 
the last section, a summary and an outlook are given. 



II. ALGORITHM 

The GS problem for general spin glasses is NP-hard2i. 
This means that only algorithms are known, where 
the running time in the worst case increases exponen- 
tially with the system size 2 ^. Hence, only small sys- 
tems can be studied in the general case. For the spe- 
cial case of two-dimensional spin glasses without an 
external field and with periodic boundary conditions 
in at most one direction, i.e. for planar graphs, ef- 
ficient polynomial algorithms 23 for the calculation of 
exact GSs are available. If one is interested in ob- 
taining the partition sum, without obtaining spin con- 
figurations, one can also treat systems with full peri- 
odic boundary conditions in polynomial time, by us- 
ing transfer-matrix approaches^SSiSSiSl, but the running 
time is again strongly increasing, limiting the investiga- 
tions to small systems. The most recent studies are based 
on matching algorithm o 2 ?i 3 ?i 31 i 32 i 33 i 34 i 3 ^i 3 ^ 3 7 , while other 
exact approaches can be found in Refs^Si 3 ^^. Results 
for systems of size 1800 x 1800 have now been obtained^! 

Here, just the basic idea of the matching algorithm will 
be explained. For the details, see Refe. HHHi The 
method works for spin glasses which are planer graphs. 
In the left part of Fig. ^ a small 2d system with open 
boundary conditions is shown. All spins are assumed to 
be "up" , hence all antiferromagnetic bonds are not sat- 
isfied. If one draws a dotted line perpendicular to all 
unsatisfied bonds, one ends up with the situation shown 
in the figure: all dotted lines start or end at frustrated 
plaquettes and each frustrated plaquette is connected to 
exactly one other frustrated plaquette. Each pair of pla- 
quettes is then said to be matched. Now, one can consider 
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FIG. 1: 2d spin glass with all spins up (left, up spins not 
shown). Straight lines are ferromagnetic, jagged lines antifer- 
romagnetic bonds. The dotted lines connect frustrated pla- 
quettes (crosses). The bonds crossed by the dotted lines are 
unsatisfied. In the right part the GS with three spins pointing 
down (all other up), corresponding to a minimum number of 
unsatisfied bonds, is shown 



the frustrated plaquettes as the vertices and all possible 
pairs of connections as the edges of a (dual) graph. The 
dotted lines are selected from the edges connecting the 
vertices and called a perfect matching, since all plaquettes 
are matched. One can assign the edges in the dual graph 
weights, which are equal to the sum of the absolute val- 
ues of the bonds crossed by the dotted lines. The weight 
A of the matching is defined as the sum of the weights of 
the edges contained in the matching. As we have seen, A 
measures the broken bonds, hence, the energy of the con- 
figuration is given by E = —J2{i j) I + 2 A. Note that 
this holds for any configuration of the spins, since a corre- 
sponding matching always exists. Obtaining a GS means 
minimizing the total weight of the broken bonds (see 
right panel of Fig.Q, so one is looking for a minimum- 
weight perfect matching. This problem is solvable in poly- 
nomial time. The algorithms for minimum-weight per- 
fect matchingsi 4 ^ are among the most complicated algo- 
rithms for polynomial problems. Fortunately the LEDA 
library offers a very efficient implementation 4 ^, except 
that it consumes a lot of memory, which limits the size 
of the systems to about N = 500 2 on a typical 500 MB 
workstation. 

For system with periodic boundary conditions in both 
directions, the matching approach is not feasible. Hence, 
we restrict the study to systems with free or with half 
periodic-half free boundary conditions. 

To study the behavior of two-dimensional spin glasses, 
here not only GSs, but also low-lying excited states are 
generated. For this purpose GS algorithms can be used 
as well. The general procedure consists of these three 
steps: 

1. Calculate the GS {sf } } of a given realization and 
the GS energy Eq. 

2. Modify some of the couplings so that the GS is 
changed. 

3. Calculate the GS {S^} of the modified system, 
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which is usually a low-lying excited state of the 
original realization. The energy of {Sl m) } calcu- 
lated with the original (unmodified) bonds is de- 
noted by Eg 1 . 

In this work four different ways of generating excita- 
tions are considered. The technical details are in the 
corresponding sections. 

• Single-bond droplets: the most straightforward im- 
plementation, where only one bond is changed with 
respect to the initial coupling. 

• Cross droplets: This mimics the generation of 
droplets induced by flipping one (central) spin. It 
involves changing O(L) bonds, iterating over 2L — 2 
changes, and picking the minimum-energy excita- 
tion among them. 

• First excitations: The lowest excitation above the 
GS is calculated. This involves the change of one 
bond (or O(L) bonds if the boundary spins are 
fixed), iterating over 0(L 2 ) changes, and picking 
the minimum-energy excitation among them. 

• Fixed-volume droplets: The matching algorithm 
does not allow for fixing the size. Hence, the size 
constraint is imposed here, when selecting among 
different cross droplets of each realization. 



The main quantities analyzed are the energy AE = 
En, the volume n and the surface S of the 



droplets/excitations. 



III. SINGLE-BOND DROPLETS 

In Refill, droplets were generated by flipping one cen- 
tral spin with respect to the GS. Since matching algo- 
rithms do not allow the inclusion of external fields, this 
approach cannot be directly applied. The most natural 
and simple idea is to flip instead a central bond with re- 
spect to the GS. This is performed by introducing a hard 
bond, i.e. changing the value of one bond {io,jo) such 
that \Ji j 1 is so large \Ji j a | = ./big, that it will be always 
satisfied in a GS, see Fig. [21 In this case the hard bond 
(*0)Jo) will be inverted, i.e. the spins Si and Sj will 
be forced to take a relative orientation opposite to the 
GS, hence Ji j = —Jhi g Si Sj Q . The GS of the modified 
system now generates a minimum-energy droplet, called 
a single-bond droplet here, with the constraint that the 
surface of the droplet runs through (io,jo). 

We have studied single-bond droplets for system sizes 
L = 6 to L = 240 while averaging over more than 3000 
samples for each size. In Fig.[3]the average droplet energy 
(AE) is shown as a function of the system size L. A 
double-logarithmic plot (not shown) shows a curvature, 
which is not compatible with a convergence of the droplet 
energy to zero, but compatible with a convergence to a 
non-zero energy. This seems plausible, because the choice 




FIG. 2: Method used to generate the single-bond droplets. Af- 
ter calculating the GS, one inverted hard bond (see triangle) 
generates an excitation (dark inner area). 



that a certain bond, i.e. the inverted hard bond in the 
modified system, has to belong to the minimum domain 
wall, imposes a penalty on the energy of the droplet. This 
energy should be of the order of the coupling constant 
J = 1 . In the spirit of Ref . ITfil where a crossover of the 
cross droplet energy according to AL e + BL^^ has been 
observed, we have performed a fit to a function AE(L) = 
AE r + AL e + BL- U with fixing 6 = -0.29, resulting in 
AE r = 1.79(6) and w = -1.8(2), with a fair quality 
of the frfcsl, Q = 0.14. Interestingly, this is a similar 
correction to scaling exponent ix> as found when iterating 
the generation of excitations over all single bonds and 
selecting the droplet from all excitations containing the 
central spini&. One might be tempted to conclude that 
these types of droplet shows a similar behavior to the 
cross droplets studied in Ref. fill (see also the inset of 
Fig. EJ. But the fit is not very sensitive to the choice 
of 9. If one allows also the variable 8 to adjust in the 
fit, 6 = -1.01(9), AE r = 1.96(2) and lu = -6(4) results 
in Q = 0.83. Hence, this type of droplet seems to be 
very different from the cross droplets studied in Ref. Il6l 
(which are similar to the droplets generated by flipping 
a single spin within a Monte-Carlo minimization^). 

This difference is visible also in the behavior of the 
single-bond droplet volume and surface, which are shown 
in Fig. 0] We have fitted to power-law functions L Dn 
resp. L ds resulting in exponents D n = 1.39(1) and ds = 

0. 769(5), hence D n /ds — 1.8. This means this type of 
droplet grows much slower than the size of the system, 

1. e. it is very different from growth as L 2 which was found 
for cross droplets^. 

Basically the above results mean that single-bond 
droplets are not the droplets which will dominate the 
low-temperature thermodynamics of the spin glass, due 
to their high energy cost as mentioned above. Therefore, 
other procedures have to be applied to obtain more phys- 
ically relevant droplets when using a matching algorithm. 
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FIG. 3: Average energy of single-bond droplets. The line 
shows a fit to the function AE(L) = AE r + AL 6 + BL~" 
with 9 = -0.29 resulting in AE r = 1.79(6) and u = -1.8(2). 
The inset shows the same data with the limiting energy AE r 
subtracted in a double logarithmic plot. 
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FIG. 4: Average volume n and average surface 5* for the 
single-bond droplets as a function of the system size L. The 
straight lines represent fits to power-laws with exponents 
D n = 1.39(1) resp. d s = 0.769(5) 



However, the single-bond droplets are not without their 
uses, as they are the first step in obtaining the lowest 
energy excitations of the system, (see section V). 



IV. CROSS DROPLETS 

The basic idea of the cross droplets is to mimic the 
approach, which was used by Kawashima and Aoki to 
obtain droplets using a Monte Carlo simulationii. They 
first calculated the GS. Then they recalculated the GS 
with the constraints that the spins on the boundary keep 
their GS orientations while a central spin is flipped. Us- 
ing this approach, small system up to L — 50 could be 
studied, and a scaling AE ~ L e of the droplet energy 
with 9' = -0.45(1) was found. 




FIG. 5: Method used to generate the cross droplets. After 
calculating the GS, several hard bonds (thick lines) are intro- 
duced, one hard bond is inverted (see triangle), leading to the 
appearance of an excitation (dark inner area). 

Our generation of the droplets works in the following 
way. After obtaining the GS several hard bonds are in- 
troduced, see Fig. If the subsystem of hard bonds 
does not exhibit frustration, no hard bond will be bro- 
ken, when a new GS is calculated. First, all boundary 
spins are fixed relative to each other by introducing hard 
bonds around the border, i.e. the bonds between pairs 
of boundary spins are replaced. The sign of the hard 
bonds is chosen such that they are compatible with the 
GS orientations of the adjacent bonds. This keeps the 
spins on the boundary in their GS orientations. Second, 
a line of hard bonds is created which runs from the mid- 
dle of (say) the left border to a pre-chosen center spin, 
again fixing the pair's spins in their relative GS orienta- 
tions. Next, the sign of exactly one hard bond on this 
line is inverted. Finally, a GS of the modified realization 
is calculated. With respect to the original GS, the result 
is a minimum energy excitation fulfilling the constraints 
that it contains the center spin, does not run beyond the 
boundary and that it has a surface which runs through 
the hard bond which has been inverted. The energy of 
the excitation is defined as the energy of the resulting 
configuration calculated using the original bond config- 
uration. For each realization, this procedure is iterated 
over all the bonds which are located on the line from 
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the boundary to the center, when in each case exactly 
one hard bond is inverted. Furthermore, this procedure 
is iterated over all four lines of bonds running from the 
left, right, top and bottom boundary to the center spin. 
Among all generated 2L — 2 excitations, the one exhibit- 
ing the lowest energy is selected as the cross droplet. 

In a previous worklfi, we have found that the cross 
droplets are almost equivalent to those of Kawashima and 
Aoki. Indeed for the system sizes which can be studied, 
the difference is much smaller than the statistical error 
bars. It turned out that cross droplets exhibit for small 
sizes a power-law behavior with an effective exponent 
6' = —0.47. But for large sizes a crossover to a power- 
law behavior is found, governed by the exponent 8 = 
—0.29, which is the same value as found in domain-wall 
studies 2 ! 4 ! 6 ! 10 . 11 ! 12 . 

Now we want to compare our results with recent re- 
sults on fixed-volume droplets obtained by Berthier and 
Young 17 . They have studied the energy dependence of 
droplets for systems with L = 64 when fixing the droplet 
size to an interval [n,n + 8(ri)]. They considered two 
cases, S(n) = and 8(n) = n — 1. For each realization, 
they optimized the energy of a droplet using a parallel- 
tempering Monte Carlo simulation. The optimization 
was done while observing the constraints that the cen- 
ter spin is included, that the size is within the given size 
range and that the droplet remains connected. For both 
cases of S(n), they observed for larger droplets a scal- 
ing of the droplet energy which is compatible with the 
DS scenario, i.e. a AP ~ n e */ d f behavior with 8* — 
around -0.23 and df = 2. For the case of fixed droplet size 
(5(n) = 0), this behavior was found also for the smallest 
droplet sizes. For the fluctuating sizes (8(n) = n — f), 
a crossover from the small-droplet behavior with a more 
negative effective exponent to the large-droplet behavior 
was found, similar to the crossover observed before for 
the cross droplets 1 ^. 

Furthermore, they have studied the distributions 
P(AP) of droplet energies. In both cases they have found 
that the distribution exhibits a scaling behavior accord- 
ing to 

For the case S(n) = they found that Pl(0) is zero (or 
very) close, while for 5(n) = n a convergence to a non- 
zero value was obtained. 

This difference between their two droplet types mo- 
tivated us to study P(AP) also for cross droplets. In 
Fig. [S] the distributions of the droplet energies for this 
droplet type is shown. We have studied system sizes up 
to L = 160 and averaged over 5000 independent realisa- 
tions in all cases. One observes that P(0) is different from 
zero and growing with system size. The typical droplet 
energy, i.e. the maximum of the distribution, is at a value 
different from zero in the thermodynamic limit. This 
can be seen from Fig. {7\ where the rescaled distributions 
(AP)P(AP) of the droplet energies for the cross droplet 



type are shown. The scaling behavior of (AE) shows 
finite-size corrections, as discussed in Ref. In that 
work, for large systems, (AE) ~ L e with 8 ~ —0.28 was 
found. Interestingly, the scaling assumption @ seems 
not to work well near AE = 0, where a systematic in- 
crease of the data points with system size is found, in- 
stead of a data collapse. To investigate this effect, we 
have evaluated P(0) as a function of system size, see in- 
set of Fig. [3 Here again a crossover can be observed. 
For larger sizes, we find P(0) ~ L e °, with 8 = 0.45(1). 
Hence the scaling assumption P^(AP) ~ -jjP(^y-) may 
be wrong. But given the system sizes studied here, it 
cannot be excluded that for even larger sizes, P(0) shows 
another crossover to the behavior found for the mean 
value (AP) already at smaller systems. 

To conclude, the scaling behavior of the distribution 
of droplet energies for the cross droplets is similar to the 
S(n) = (n — 1) droplets of Berthier and Young, in the 
sense that P(0) is finite and growing. This similarity has 
been observed already, when studying simply the mean 
value AP. On the other hand, the scaling assumption 
(PJ seems not to work near AP = in contrast to the 
results of Ref. ^3- Thus, we have another example of the 
scaling behavior of droplets seemingly depending on the 
recipe used to generate them, at least at the system sizes 
which can be studied at present. 




1 2 3 4 5 

AE 



FIG. 6: Distribution of droplet energies for cross droplets for 
system sizes L = 14, 20, 40, 80, 160. Error bars are of the order 
of the symbol size (see L = 14 datapoints). Lines are guides 
for the eyes only. 



V. FIRST EXCITATIONS 

Next, we consider the first (i.e. the lowest) excitation 
of the system above the ground state. They have been 
studied by Picco, Ritort and Sales 19,20 ' who generated the 
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FIG. 7: Rescaled distribution of droplet energies for cross 
droplets for system sizes L — 14, 20, 40, 80, 160. The inset 
shows the probability P(0) to get a zero-energy droplet as a 
function of system size L. For larger sizes, P(Q) scales as L e ° 
with 6» = 0.45(1). 



lowest excitation for sizes up to L — 16 using a transfer 
matrix method. They have measured the exponent 9\ 
describing the finite-size scaling AE ~ L e± of the ex- 
citation energy and the exponent A; characterizing the 
size distribution of the excitation volumes n according 
to P(n) ~ n~( 1+x '\ Furthermore, they have derived a 
scaling relation 9 — 9\ +dXi (d = 2 here), connecting the 
exponents to the usual droplet exponent 9. Picco et al. 
argue 9\ = —2 (while their numerical results are close to 
9\ = —1.7(1), but they explain this discrepancy as due 
to the small sizes they could study) and they have found 
A* = 0.77(1), leading to 6 = -0.46(2). 

As we have already seen for the single-spin induced 
(i.e. here cross) droplets, large sizes may be needed to 
see the correct scaling behavior, hence L = 16 may be 
much too small. To test their results for larger sizes up to 
L — 64 we have also calculated first excitations by using 
the matching algorithm. They are obtained by iterat- 
ing one hard bond over all bonds and then selecting for 
each realization the minimum excitation among all 0(L 2 ) 
droplets. We have considered free and fixed boundary 
conditions. In the latter case, again a closed loop of hard 
bonds winds around the system. We averaged over a 
number of realizations ranging from a few hundred for 
the largest sizes to 10 4 for L < 12. In Fig.[S]the average 
energy of the lowest excitations for both types of bound- 
ary conditions are shown. For free boundary conditions 
excitations leaving the boundary spins unchanged can be 
realized. Thus, for each system, the energy of the low- 
est excitation with fixed boundary conditions is an upper 
bound for the energy of the excitation with free boundary 
conditions. Furthermore, in the thermodynamic limit, 



the boundary conditions are expected to be irrelevant 
for the energies of a droplet, hence both cases should 
agree for L — > oo. Indeed in Fig. |H| the average energy 
for the fixed boundary conditions is always above the 
average energy for droplets with free boundaries. Fur- 
thermore, both curves approach each other with growing 
system size. For small system sizes, the first excitations 
with free boundary conditions can take advantage of the 
system boundary, i.e. many of the droplets will run up 
to the boundary, where creating a domain wall does not 
cost energy. Hence large corrections to scaling can be 
expected, as is visible in the figure. In Ref. only small 
systems L < 12 were considered. Hence, the corrections 
to scaling were hardly visible and a i -17 behavior was 
observed. 

Droplets for the case with fixed boundaries show less 
corrections to scaling, because the system cannot create 
domain walls of low energy cost. For L > 10 a fit to an 
algebraic decrease of the excitation energy of the form L b 
yields b — —2.27(2). Since the lowest excitations for both 
boundary conditions have to agree for L — > oo and be- 
cause the droplets with fixed boundaries exhibit almost 
no corrections to scaling, it can be expected that also the 
lowest excitations with free boundaries show the L b scal- 
ing for large sizes. This can be represented by a function 
with al 6 /(1 + eL b ~ b ) scaling behavior, see Fig. [SJ But 
the system sizes accessible to us, although much larger 
than those of Ref. IT^. are still too small to give a definite 
answer for this type of droplet. 




FIG. 8: Average excitation energy AE as a function of system 
size for lowest excitations with free and with fixed boundary 
conditions. The lines denote fits to functions of the form 
A£ flxcd (L) = aL Bl with a = 55(4), 0i = -2.27(2) (fitted for 
L > 10) and A£ frcc (L) = cL b ' /(l + eL 6 '" 91 ) with c = 2.8(2), 
b' = -1.27(5), e = 0.041(3). 

The fact that the droplets with free boundaries have 
larger finite-size corrections can be seen as well when 
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studying the distribution of the droplet energies. In Fig. 
ED the rescaled distribution P(AE)L ei of the droplets is 
displayed. For small systems a good data collapse is ob- 
tained for B\ = —1.4, corresponding to the slope of the 
mean droplet energy in Fig. [5] For larger sizes, the best 
data collapse is obtained for 9\ = —1.8. 
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FIG. 9: Rescaled distributions P(AE)L Bl as a function of 
the rescaled excitation energies AE/L 1 for the lowest exci- 
tations for free boundaries. For large system sizes 0\ = —1.8 
was used, while for small systems (see inset) 9i = —1.4 was 
applied. As an example, the error bars for L — 32 are shown, 
which are large in the region of large, hence rare, energy val- 
ues. 

For the case of fixed boundaries, corresponding to the 
Fig. El the best data collapse P(AE)L ei versus AE/L 91 
is obtained for B\ = —2.3, see Fig. 1101 

The different values for 0\ for the L 01 scaling can be 
understood as follows. The lowest excitation is generated 
by iterating one hard bond over all bonds of a system. 
Then the minimum-energy droplet is taken among all this 
0(L 2 ) excitations. To the lowest order, this procedure 
generates droplets, where the energy scales like 



AE~ 1/(L 2 P*(0)), 



(4) 



here P*(0) is the probability to obtain a (single-bond) 
excitation having zero energy. Hence, to understand the 
scaling of AE, one should study the scaling of P*(0). 
In Fig. ^2 the distribution of energies of all single- hard- 
bond excitations (hxed boundary conditions) is displayed 
along with the finite-size behavior of P*(0). For large 
enough sizes P*(0) scales about as L 3 leading to AE ~ 
L~ 2 - 3 as found in the data. Hence, the scaling P*(0) for 
the single-bond excitations is compatible with a P*(0) ~ 
1/L behavior with 6 rj —0.3, yet again the value found 
from domain wall studies. 

Note that for small sizes, P*(0) grows more rapidly 
than at larger sizes, while from the scaling behavior of 
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FIG. 10: Rescaled distributions P(AE)L Bl as a function of 
the rescaled excitation energies AE/L 1 for the lowest excita- 
tions for fixed boundaries. 6\ = —2.3 is used. As an example, 
the error bars for L = 20 are shown. 



(AE), the opposite would be expected. This deviation 
is presumably due to the fact that for small sizes the 
behavior of the energy distribution a little bit away from 
AE = contributes as well. In this region it was assumed 
when deriving Eq. that the distribution is constant, 
so (P*(AP*) ~ P*(0)) for small energies, which is, as 
visible in Fig. ^] not the case. It is only for larger sizes 
that the scaling behavior at AE — is relevant, leading 
to the expected result. 

In Fig. E| the corresponding distribution of ener- 
gies for the excitations with free boundary conditions 
is shown. As one can see, excitations having zero en- 
ergy are much more likely than in the fixed boundary 
case, because they can take advantage of letting the do- 
main wall end at the boundary. For small sizes P*(0) 
is slightly decreasing, corresponding to the moderate de- 
crease of (AE) seen in Fig. [H] while for moderate sizes, 
P*(0) is constant as a function of system-size, indeed 
compatible with 6\ = —2, as argued by Picco et al. Note 
that here the scaling of (AE) is also for small sizes com- 
patible with the scaling of P*(0), because here P*(AP*) 
is much more constant near AE* = and P* (0) is much 
larger. 

For larger system sizes, when the system boundary 
is far away, the behavior of P*(0) for the two cases 
of boundary conditions should agree again. Unfortu- 
natly, although the matching algorithm is a very powerful 
tool, these sizes are out of reach with current technology. 
Hence it is hard to say, whether P*(0) ~ L 0,3 or P*(0) 
constant is the true limiting behavior 

Finally, to complete our comparison with the results 
of Ref. |2(J we show in Fig. ^| the distributions of exci- 
tation volumes for the lowest excitations for both types 



8 




FIG. 11: Distribution P*(AE*) of single- bond excitations en- 
ergies with fixed boundary condition when iterated over all 
bonds of each realization, for system sizes L = 4, 10, 20, 48. 
The inset shows the finite-size behavior of P*(0). The line in 
the inset represents a L 0,3 behavior. 
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FIG. 12: Distribution P*(AE*) of single-hard-bond exci- 
tations energies with free boundary conditions when iter- 
ated over all bonds of each realization, for system sizes 
L = 4, 16, 32, 64. The inset shows the finite-size behavior 
of P*(0). 

of boundary conditions. It exhibits a P{n) ~ n~^ l+ l ' 
behavio»2£, except for the largest sizes, where the excita- 
tions start to interact with the system boundary. We 
have fitted power-laws to our data, ignoring the data 
points for the largest sizes and obtained A; = 0.67(4) 
for free be and A; = 1.04(5) for fixed be. The former 
value is compatible with the previous result A; ~ 0.7 



FIG. 13: Distributions P(n) of excitation volumes for lowest 
excitations with free boundary conditions for different system 
sizes L = 16, 24, 32, 48, 64. Apart from the largest volumes, 
where the excitations start to touch the boundary, the data 
is very well described by a power-law behavior n~ 1 ~ x ' with 
A; = 0.67(4) from a fit to the L — 48 data in the range 
n € [1, 1000]. The inset shows the same distribution for the 
case of fixed boundary conditions. The line represents in this 
case a power law with A; = 1.04(5)). 

obtained from direct finite-size scaling^, while the final 
value quoted in this reference is X eS = 0.77 obtained from 
an aspect-ratio analysis, i.e. using non-square systems. 

Anyway, it seems strange that in both cases no large 
size-dependence of A; is visible. Since one would again 
expect that in the thermodynamic limit results obtained 
from both types of boundary conditions should agree, one 
would expect that the behavior of one of the distributions 
should converge to the behavior of the other. 

To summarize, for the behavior of the mean value 
(A£J), the free-boundary case exhibits a strong size de- 
pendence, while the fixed-boundary case exhibits almost 
none. Hence in this case one would expect that for large 
sizes the free-boundary case seems converged to the fixed- 
boundary results. When studying P*(0), both cases show 
a crossover, hence it seems possible that one or both 
exhibit another crossover, leading to the same thermo- 
dynamic limit. Finally, for the distribution of sizes, it 
is hard to imagine that either shows a strong crossover. 
There is a possibility that the boundary conditions may 
indeed matter, and that a convergence to the same be- 
havior may not take place. 



VI. FIXED- VOLUME DROPLETS 

To compare with the results of Berthier and Young 17 
we have implemented a study similar to theirs for the 
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FIG. 14: Average droplet energy for cross droplets (L = 160), 
when the selection of the minimum droplet is restricted to a 
excitation size interval [n, n + S(n)]. Two cases S(n) = 9 and 
8(n) =ri—l are shown, in the spirit of the droplet generation 
mechanism in Ref. 



size dependence for the cross excitations. For the largest 
system size L = 160, for each realization and a given size 
window [n, n + 5(n)], the minimum-energy droplet was 
selected only among the cross excitations which lie inside 
the size window. Here we have studied also the case 
6(n) = n — 1. Furthermore, instead of taking S(n) = 0, 
we have taken S(n) = 9 to improve the statistics. 

In Fig. 1141 the resulting average droplet energies are 
shown. However, the behavior we find here is totally 
different from that reported in Ref. Il7l For the window 
of fixed length, the average droplet energy even goes up 
with n, while for the case 5(n) = n— 1, the droplet energy 
decreases, but for small sizes with a slower n~ 53 power- 
law decrease than the n -0 ' 22 decrease of Ref. [lTj For 
larger droplet sizes, a crossover can again be observed, 
and the data is compatible with a n~°' n decrease, which 
compares well with the M " 0115 behavior found in Refill 
For very large sizes n, the excitations touch the boundary 
of the systems, where the spins are fixed. This leads to 
a strong increase of the droplet energy. 

To understand our result for small and medium 
droplets, one should first note that there is an impor- 
tant difference between the droplets beind studied here 
and in Ref. ^3- There, first the size was fixed, then the 
optimization was performed while obeying the size con- 
straint. Here, first for each given inverted hard bond, a 
first optimization is performed, with a restricted search 
space. Then the sizes of the resulting excitations are mea- 
sured and finally the minimum-energy droplet is selected 
among the excitations having the wanted sizes, i.e. a sec- 
ond optimization is performed. Note that the matching 
algorithm does not allow for fixing the size of a droplet 



in advance, so one cannot circumvent the two-stage op- 
timization process. 

We assume now that the behavior of AE can be under- 
stood in the following way. The distribution of excitation 
energies P n (AE*), obtained after the first minimization 
procedure, is taken as input. For each size, the final 
droplet is selected as minimum among a given number of 
excitations which we denote as M(n). Hence a behavior 
AE ~ l/(M(n)P n (0)) can be expected, as in the case of 
the first excitation discussed in the last section. We shall 
now test this assumption using the available data for the 
case S(n) = n — 1. 
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FIG. 15: Distribution of energies for all cross excitations (L = 
160), restricted to a size interval [n,n + 5(n)], here 5(n) = 
n — 1. The inset shows the behavior of P(0) as a function of 
size. The line in the inset denotes a power law n -0,49 

In Fig. ^] the distributions P n (AE*) of excitation 
energies are shown for different sizes n. Also the de- 
pendence of P„(0) on n is shown in an inset. Due to 
numerical problems in determining P„(0), especially for 
small excitations, the fluctuations are quite large. If one 
assumes a power law decrease, the data is roughly com- 
patible with a n~ - 49 behavior for small n. But, in fact, 
no region with a clear power-law behavior can be easily 
identified here. 

M(n) can be obtained from the distribution P{n) 
of excitation sizes, see the inset of Fig. ^1 via 
M(n) r~ f ™ +S ^ P(n)dn. In the main part of Fig. ^3 
c/(M(n)P n (0)), with c = 2700 chosen by hand, is com- 
pared to AE. The agreement is fair, and the large fluc- 
tuations originate in the numerical problems when deter- 
mining P n (0). One can circumvent this fluctuations for 
the small-size data, by using scaling arguments. From 
a fit to the small size range n £ [10, 100] one obtains 
P(n) ~ n-°- 46 ( 2 ). Since the selection is done in a range 
0(n), one has M(n) ~ n~ 0A6+1 ~ n 54 . In combina- 
tion with the above result P n (0) ~ n -0 ' 49 one obtains 



10 




10 100 1000 10000 



n 

FIG. 16: Average droplet energy for cross droplets (L = 160), 
when the selection of the minimum droplet is restricted to a 
excitation size interval [n, 2n — 1]. The numerical data from 
Fig. 1141 is repeated here ("measured") and compared to the 
function c/(M(n)P n (0)), with the scaling factor c = 2700 
("indirect"). The inset shows the distribution of sizes of the 
cross excitations. M(n) is obtained by integrating P(n) over 
intervals of size S(n). 



AE ~ l/(M(n)P n (0)) = l/(n a54 n-°- 49 ) = n" ' 05 , in 
agreement with the direct fit of the numerical data, see 
Fig. m 

For the case of fixed window size, each droplet is se- 
lected among M(n) ~ P(n) ~ n ~0A6 d r0 pi e ts. Hence a 
increase of AE is not surprising here. But the numer- 
ical problems due to a much smaller (0(1) opposed to 
0(n)) statistics when determining P n (0) are too large in 
this case, such that a direct comparison cannot be per- 
formed. 

To summarize, the two recipes to generate fixed-size 
droplets are different here and in the work of Berthier 
and Young. Hence, in general different results come out. 
But for the case of size fluctuations S(n) = (n — 1), the 
scaling behavior is the same. Since fluctuations 0(n) of 
the droplet size seem to be the most natural case, directly 
in the spirit of the droplet pictures by Fisher and Huse a , 
it seems not surprising that in this case the results are 
similar. In our case it stems from a complex interplay of 
the selection mechanism and the behavior of the distri- 
bution P(n) of excitation sizes as well as the distribution 
P n (AE*) of excitation energies for small energies. No 
direct relation to the scaling exponent 9 is obvious. Fur- 
thermore, unfortunately, this still does not allow us to 
understand why for the S(n) = 0(1) case Berthier and 
Young found the same scaling behavior, but without the 
crossover usually seen at small to medium system sizes. 



VII. SUMMARY AND OUTLOOK 



We have studied droplets and first excitations for two- 
dimensional spin glasses with Gaussian interactions. Us- 
ing highly sophisticated matching algorithms, large sizes 
can be studied, but one is limited in the flexibility of the 
recipes for droplet generation. 

We studied four different types of droplets/excitations 
and observed in general that the results depend strongly 
on the way droplets are generated. The single-bond 
droplets were very simple to generate. They do not com- 
pare with typical droplets, which are assumed to be ob- 
tained e.g. by the method of Kawashima and Aoki, be- 
cause in the thermodynamic limit their energy converges 
to a non-zero value. The cross droplets, already intro- 
duced in Ref. 2 are almost equivalent to the single- 
spin droplets, i.e. can be considered as typical, and 
they show many properties expected from scaling the- 
ory. Much larger sizes than with Monte Carlo methods 
can be studied, hence the crossover to the L 6 {9 = —0.29) 
behavior observed. Also the distribution of the droplet 
energies is almost as expected, except the fact the the 
scaling assumption does not work very close to zero en- 
ergy. When studying first excitations, again much larger 
sizes than in previous works are possible when using the 
matching approach. Different results are obtained for 
free and fixed boundary conditions. Unfortunately, in 
this case the range of accessible sizes is even for the pow- 
erful algorithm used still too small to decide whether in 
the thermodynamic limit both cases agree and what the 
final scaling behavior will be. Finally, the study of fixed- 
volume droplets shows the limitations of the matching 
approach. It is not possible to use the matching algo- 
rithm to generate droplets that are really equivalent to 
the droplets generated by Berthier and Young. Never- 
theless, it is possible to understand the scaling behavior 
of the droplets generated here in terms of the distribu- 
tions of the excitations and their sizes. Furthermore, for 
the case of fluctuating volumes, which is believed to rep- 
resent the typical behavior in the spirit of the DS theory, 
even the same scaling behavior as obtained by Berthier 
and Young is found. 

For the single-bond and the cross droplets, their be- 
havior seems to be clear. For the first excitations, since 
the matching algorithm runs in polynomial time, and be- 
cause the generation of first excitations requires 0(L 2 ) 
calls, the arrival of more powerful computers will resolve 
sooner or later the question concerning the role of bound- 
ary conditions. It would be desirable to tackle also the 
last open point raised by this paper: is it possible to ex- 
tend or modify the matching technique, such that real 
fixed-size droplets can be studied? This would allow to 
understand their behavior better as well. 
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